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OPERATING  CHARACTERISTICS  FOR  WEIGHTED 
ENERGY  DETECTOR  WITH  GAUSSIAN  SIGNALS 

INTRODUCTION 

The  operating  characteristics  of  an  equi-weighted  energy 
detector  for  Gaussian  signals  in  noise,  in  terms  of  false  alarm 
and  detection  probabilities,  can  be  characterized  mathematically 
by  a  partial  exponential  expansion,  and  have  previously  been 
numerically  evaluated  for  arbitrary  numbers  of  samples  and 
signal-to-noise  ratios  11;  (7)  -  (8)  and  figures  2-6}. 

However,  when  the  weights  employed  in  the  energy  detector  are 
unequal,  or  if  the  signal  and  noise  powers  on  each  sample  are 
unequal,  these  results  do  not  apply  and  can  be  misleading, 
especially  when  the  number  of  samples  summed  is  not  large,  what 
is  needed,  in  this  case  of  arbitrary  numbers  of  samples  and 
unequal  weights  or  powers,  is  an  exact  approach  in  terms  of  the 
characteristic  function  of  the  decision  variable;  this  latter 
function  is  frequently  available  in  closed  form  and  can  be 
employed  in  the  fast  efficient  procedure  presented  in  12)  and 
utilized  in  13,4,5}  for  direct  accurate  evaluation  of  the 
exceedance  distribution  function. 

At  the  same  time,  it  would  be  very  useful  to  have  accurate 
approximations  for  the  receiver  operating  characteristics,  which 
apply  over  the  full  range  of  applicable  false  alarm  and  detection 
probabilities,  yet  are  easily  computed  in  terms  of  readily 
available  functions,  or  circumvent  some  of  the  more  difficult 
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numerical  procedures  required  ir.  the  exact  approach.  Here,  we 
will  consider  four  such  approximations,  namely  Gaussian,  chi- 
square,  constant  plus  chi-square,  and  generalized  noncentral 
chi-square,  and  demonstrate  the  range  of  applicability  of  each. 

Thus,  our  goals  here  are  two-fold 

(1)  determination  of  exact  operating  characteristics  of 
arbitrary  weighted  energy  detectors  along  with  working  programs, 
thereby  allowing  for  investigation  of  other  similar  cases  of 
interest  to  the  user;  and 

(2)  construction  of  accurate  simple  approximations  to  the 
operating  characteristics,  which  can  be  extended  to  related 
difficult  problems  and/or  circumvent  complicated  numerical 
procedures . 

As  a  by-product,  the  inadequacy  of  some  extant  approximations 
will  be  delineated  quantitatively;  in  particular,  the  generally 
optimistic  results  predicted  by  the  Gaussian  approximation  will 
be  shown  to  prevail  even  when  the  number  of  independent  samples 
involved  in  the  energy  detector  is  very  large. 
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CHARACTERISTIC  FUNCTION 

We  presume  that  we  have  M  channels  (or  samples)  containing 
either  noise-alone  or  signal-plus-noise,  and  that  the  random 
variables  in  each  channel  are  statistically  independent  of  each 
other.  Specifically,  for  our  interest,  the  output  envelopes, 

{effl}  for  1  <  m  <  M,  of  M  disjoint  narrowband  filters  are 

•* 

subjected  to  weighted  square-law  summation  for  purposes  of 
threshold  comparison  and  a  statement  about  signal  presence  or 
absence  on  that  particular  observation  of  M  outputs.  The 
decision  variable  in  this  case  is 

M  M 

E2  \  ' 

%  em  ’  *»  '  1 

m»l  m«l 


where  weights  {wm}  are  all  positive  but  otherwise  arbitrary,  and 
the  M  squared-envelope  outputs  {zm)  are  statistically  independent 
and  identically  distributed.  An  example  is  afforded  by  a  finite¬ 
time  exponential  summer  where  wffl  «  A  rm_*  ,  r  <  1,  1  <  m  <  M. 

Without  loss  of  generality,  the  sum  of  the  weights  is  set 
equal  to  unity, 


that  is,  A 


1-r 


(2) 


Then,  the  mean  of  random  variable  x  in  (1)  is  equal  to  the  mean 
of  each  random  variable  z  .  because  all  the  {z  1  are  identically 
distributed.  (If  there  are  scaling  differences  in  the  variables 
(zm),  these  factors  can  be  absorbed  in  modified  scalings  {w^} , 
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without  loss  of  generality.)  Under  these  assumptions,  it  should 
be  observed  that  the  performance  of  the  weighted  energy  detector 
in  (1)  is  independent  of  the  ordering  of  the  weights;  thus,  the 
weights  can  be  arranged  in  any  order  without  affecting  the 
detection  capability.  Also,  the  absolute  level  of  the  {wm) 
cannot  affect  the  operating  characteristics  of  detector  (1). 


STATISTICS  OF  Zm 

m 

For  Gaussian  signals  and  noises  present  at  the  inputs  to  the 
M  narrowband  filters  in  (1),  the  probability  density  function  of 
each  filter  output  envelope-squared  random  variable  zm  is 

P2<u)  -  j*  exp  for  u  >  0  ,  (3) 

where  parameter 


a 


1 

+  R 


for  noise-alone 

for  signal-plus-noise 


(4) 


Here,  we  have  normalized  according  to  the  noise  power;  that  is, 

the  mean  of  random  variable  zm  is  set  equal  to  1  for  noise-alone. 

This  presumption  is  equivalent  to  having  knowledge  of  the  average 

noise  level  in  the  absence  of  signal  and  can  be  accomplished  in 

practice  by  monitoring  the  filter  outputs  over  a  sufficiently 

long  past  interval  of  time.  Also,  R  is  the  signal-to-noise  power 

ratio  per  sample  at  the  output  of  each  filter. 

The  characteristic  function  of  each  random  variable  z  in  (1) 

m 

is  given  by  expectation  (ensemble  average) 
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fz(£.)  «*  E  { exp(  i  fcz ) }  -  J  du  exp(i£u)  pz(u)  «  -1  /  ^  ,  (5) 

where  we  used  (3).  The  cumulants  {X„(k)}  of  s_  are  immediately 

Z  111 

available  from  (5)  as 

ffclYjT  X2(k)  -  ak  for  k  >  1  .  (6) 

Actually,  these  are  scaled  cumulants,  by  the  factor  l/(k-l)l; 
they  are  more  convenient  and  will  ba  employed  henceforth. 

CHARACTERISTIC  FUNCTION  OF  OUTPUT  x 

The  characteristic  function  of  summation  random  variable  x  in 
(1)  is  given  by  expectation 

M  I"  M  l”1 

fx(£)  -  E{exp(Ux)}  -  y~j  f2(waO  -  j~J(l  -  i  £>ma  j  ,  (  7 ) 

m»l  lm*l 

where  we  used  the  4 ~ dependence  of  the  {zBJ  and  relation  (5).  The 
(scaled)  cumula  .  x  are  available  from  (7)  according  to 

K 

(TTTrr  Vk>  -  »k  * ak  wk  £or  k  * 1  ■  "» 

m«l 

In  particular,  the  mean  and  variance  of  x  are,  upon  use  of  (2), 

ux  -  Xx(l)  -  a  wx  -  a  ,  *  xx(2)  ^  a2  W2  .  (9) 

The  desired  closed  form  for  the  characteristic  function  of  x  is 

given  by  (7),  where  the  signal-to-noise  ratio  parameter  a  is 
given  by  (4).  Result  (7)  applies  for  arbitrary  M,  weights  {wffl} , 
and  per-sample  signal-to-noisv  ratio  R. 
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SOME  RELATED  RESULTS 

Characteristic  functions  of  the  form  of  (7)  occur  in  numerous 
problems.  For  example,  the  stability  of  a  spectral  estimation 
technique  employing  overlapped  FFT  processing  of  windowed  data 
encountered  this  form  {6;  (35)  and  (15)],  where  weights  were 

proportional  to  the  eigenvalues  (X^J  of  a  normalized  covariance 
function.  Another  example  is  furnished  by  diversity  combination 
in  a  channel  subject  to  partially-correlated  signal  fading;  see 
(7;  ( D-l 4 ) ] ,  (8;  (24)],  and  19].  In  particular,  the  exact 
characteristic  function  in  17]  and  (8]  took  the  form 


where  (Xm]  are  the  eigenvalues  of  a  covariance  matrix.  Parameter 
D  was  the  order  of  diversity  in  [7],  but  was  equal  to  1  in  [8], 

GAUSSIAN  APPROXIMATION  TO  EXCEEDANCE  DISTRIBUTION 

For  the  general  characteristic  function  given  by  (7)  and  (4), 
a  Gaussian  approximation  to  the  probability  density  and 
exceedance  distribution  functions  is  given  in  appendix  A.  It  is 
derived  for  arbitrary  M,  weights  (wml #  and  signal-to-noise  ratio 
R.  However,  its  applicability  to  numerical  evaluation  of 
receiver  operating  characteristics,  in  the  form  of  detection 
versus  false  alarm  probabilities,  will  be  shown  to  be  rather- 
limited  in  the  next  section. 
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EXCEEDANCE  DISTRIBUTION  FOR  ALL  WEIGHTS  EQUAL 


In  this  section,  the  weights  (wm)  in  (1)  and  (2)  are  e.,ual: 


w«  —  for  1  <  m  <  M  . 
m  n 


(11) 


The  characteristic  function  in  (7)  then  becomes 


fxU>  -  <1  -  iU/M)  "  . 


(12) 


This  corresponds  to  a  multiple  of  a  chi-squared  random  variate 
with  2M  degrees  of  freedom.  The  corresponding  probability 
density  function  is 


Px(u> 


uM“X  exp( -uM/a ] 
(M-l)i  ( a/M ) M 


for  u  >  C  , 


(13) 


while  the  exceedance  distribution  function  is,  for  u  >  0, 


Qx(u)  -  J  dt  Px(t)  -  exp(-u«/a)  e^^uM/a)  *  E^^tuM/a)  .  (14 


Here,  en(x)  is  the  partial  exponential  function  (10;  (6.5.11)), 


vx) 


U  , 


(15) 


and  we  have  defined  auxiliary  function 


E  ( x )  -  exp(-x)  e  (x)  for  x  >  0 
n  n 


(16) 


If  threshold  value  T  is  used  for  comparison  with  output  x  of 
the  energy  detector  (1),  then  the  false  alarm  probability  Pp  is 
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PF  “  °x(Tf  a-1)  "  EM-X(TM)  *  U7) 

Similarly#  the  detection  probability  PD  is,  from  (14)  and  (4), 

PD  -  QX(T,  a*l+R)  -  .  (16) 

When  T  is  eliminated  between  . (17)  and  (18),  the  operating 
characteristics  (?D  versus  Pf)  can  be  plotted,  with  signal-to- 
noise  ratio  R  as  a  parameter.  Separate  plots  are  required  for 
different  values  of  M,  the  number  of  envelope-squared  samples. 


GRAPHICAL  RESULTS 

The  receiver  operating  characteristics  (ROC)  for 

«  -  1,  2,  4,  8,  16,  32,  64,  128,  256,  512,  1024  (19) 

are  plotted  in  figures  1  through  11,  on  normal  probability  paper, 
for  false  alarm  probabilities  ranging  from  IE-10  to  .1  and  for 
detection  probabilities  ranging  from  .01  to  .999.  Signal-to- 
noise  ratios  (in  decibels)  have  been  chosen,  typically,  to  cover 
Pp,PD  possibilities  from  low-quality  pair  .01,. 5  to  high-quality 
pairs  in  the  neighborhood  of  IE-10,. 99. 

Superposed  in  figure  3  (in  dashed  lines)  is  the  Gaussian 
approximation,  for  M  •  4,  to  the  exact  exceedance  distribution 
function  in  (14);  see  appendix  A.  Three  selected  values  of 
signal-to-noi se  ratio  R  are  indicated,  namely  R  ■  4,  8,  and  12 
dB.  They  are  identified  by  a  black  dot  where  they  cross  the 
exact  operating  characteristic  for  the  same  signal-to-noi se 
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ratio.  It  is  seen  that  the  Gaussian  approximation  is  virtually 
useless  at  this  low  value  of  M,  the  number  of  samples. 

This  superposition,  of  three  representative  curves  afforded 
by  the  Gaussian  approximation,  is  continued  up  through  M  *  1024 
in  figure  11.  Again,  agreement  with  the  exact  results  is 
generally  quite  poor.  Even  at  M  -  1024,  the  required  signal-to- 
noise  ratio  from  the  Gaussian  approximation  for  PF  ®  IE-10, 

PD  -  .3,  for  example,  is  in  error  by  .3  dB. 

Furthermore,  it  should  be  observed  that  the  Gaussian 
approximation  is  always  optimistic  in  the  useful  range  of  the 
operating  characteristics;  this  bias  is  misleading  in 
quantitative  performance  predictions  applied  to  practical 
detection  systems.  Additionally,  the  case  in  this  section, 
namely  equal  weights,  is  the  most  favorable  situation  for  the 
Gaussian  approximation  to  apply  in;  any  other  distribution  of 
weight  values  makes  the  effective  number  of  weights  (Me  in  (A-6) 
and  sequel)  less  than  M,  thereby  deviating  even  further  from  an 
accurate  application  of  the  central  limit  theorem.  The  message 
to  be  conveyed  here  is  that  the  performance  capability  of  energy 
detectors  for  Gaussian  signals  and  noises  should  be  based  on 
something  other  than  the  Gaussian  approximation. 


9 


Probability  of  Detection 


7R  8753 


Probability  of  False  Rlarm 

Figure  1.  ROC  for  M-l,  Equal  Weights 
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Probability  of  False  Alarm 
Figure  2.  ROC  for  M=2 ,  Equal  Weights 
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Probability  o-f  False  fllarm 
Figure  3.  ROC  -for  MM ,  Equal  Weights 
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Probability  o-f  False  Rlarm 
Figure  4.  ROC  -for  71*8,  t’qual  Weights 
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Probability  of  False  Rlarm 
Figure  6.  ROC  for  M=32}  Equal  Weights 
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Probability  of  Detection 
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Probability  of  False  Rlarm 

Figure  8.  ROC  for  M=128,  Equal  Weights 
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Probability  of  Detection 


Probability  of  False  Rlarm 

Figure  9.  ROC  for  M=256 ,  Equal  Weights 
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Probability  o-f  False  Rlarm 
Figure  10.  ROC  for  M=512,  Equal  Weights 
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Probability  of  False  filarm 

Figure  11.  ROC  for  M=1024,  Equal  Weights 
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EXCEEDANCE  DISTRIBUTION  FOR  ALL  WEIGHTS  DIFFERENT 


In  this  section,  we  confine  attention  to  the  case  where  all 
the  weights  {wmJ  are  different  from  each  other;  that  is. 


wm  #  Wj,  if  m  *  k  ;  >  0 

m  k  m 


(20) 


Then,  we  expand  the  characteristic  function  of  x  in  (7)  in  a 
partial  fraction  expansion  according  to 


fx(U 


ftfr  -  t  ' 

Ln-l  J  m 


(21) 


where  coefficients 


m  M 


for  1  <  m  <  M  , 


FT  (w»  -  wx) 


(22) 


depend  only  on  weights  {wffi}  and  not  on  signal-to-noise  ratio  R. 

The  probability  density  function  of  x  is  then  immediately 
available  from  (21)  as 


n 

pv(u)  -  )  A  B  exp(-A  u)  for  u  >  0  , 

x  /  ,  mm  m 


(23) 


where  A„  «  l/(w  a).  The  corresponding  exceedance  distribution  is 
ram  r  ^ 


Qx(u)  -  J  at  px(t) 


e*P'-V! 
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If  threshold  T  is  used  as  the  basis  of  comparison  for  output 
* 

x  of  the  weighted  energy  detector  in  (1),  the  false  alarm  and 
detection  probabilities  follow  from  (24),  respectively,  as 

PF  "  Qx(T;  a  "  '  P0  -  QX(T,-  a  -  1  +  R)  .  (25) 

As  an  example,  if  M  -  1,  then  w^  -  1,  »  1/a,  «  1,  and 

(24). yields  Qx(u)  ■  exp(-u/a)  for  u  >  0.  Then,  (25)  gives 

r. T  N  TIr  fln  Pf1 

PF  -  exp(-T)  ,  PD  «  exp [i+rJ  "  PF  *  exp  ~r~  *  {26) 

> 

For  this  special  case  of  M  *  1,  threshold  T  can  be  eliminated  and 
PD  expressed  explicitly  in  terms  of  P^  and  R. 

GRAPHICAL  RESULTS 

The  particular  case  of  unequal  weights  that  we  shall 
concentrate  on  here  is  a  set  of  exponential  weights 

wm  -  a  r31”1  for  1  <  m  <  M  ,  r  <  1  ,  (27) 

where  scale  factor  A  is  selected  for  normalization  of  the 
weights,  according  to  (2).  Of  course,  the  absolute  level  of  the 
weights  does  not  affect  the  operating  characteristics. 

In  figure  12,  the  ROC  for  M  »  4  and  r  »  .99  is  plotted,  as 
determined  from  (25)  and  (24).  Since  r  is  close  to  1  for  this 
example,  the  weights  (27)  are  all  nearly  equal,  causing  some  of 
the  coefficients  (Bm)  in  (22)  to  be  rather  large,  in  the  range  of 
+.5E6.  This  leads  to  round-off  error  in  sum  (24)  for  the 
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exceedance  distribution  function  and  the  possibility  of  useless 
numerical  results;  however,  because  M  *  4  is  a  small  number,  the 
round-off  error  does  not  yet  show  up  in  figure  12. 

When  M  is  increased  to  3  in  figure  13  and  r  is  kept  at  .99, 
coefficients  {Bm}  in  (22)  reach  values  in  the  range  of  +.7E12, 
and  round-off  error  begins  to  show  up  as  wiggly  lines  in  the 
higher  detection  probability  values  near  .999.  We  are  using  a 
computer  with  64  bits  per  word,  which  yields  approximately  15 
decimals  of  accuracy  for  the  mantissa.  Although  coefficients 
{BmJ  can  be  calculated  very  accurately  from  (22),  they  alternate 
in  sign  and  can  be  very  large.  Then  Qx  in  (24)  requires 
differencing  of  large  numbers,  with  an  attendant  possibly 
damaging  loss  of  accuracy,  especially  for  small  Pj.. 

When  M  is  increased  by  one,  to  9  in  figure  14,  and  r  is 

maintained  at  .99,  round-off  error  is  now  significant  at  the 

upper  edge  of  the  ROC,  although  useful  characteristics  are  still 

available  for  lower  values  of  PD.  The  reason  for  this  problem  is 

that  all  the  weights  are  close  to  each  other;  in  fact,  the  M-th 
1 

weight  is  r  -  .923  times  as  large  as  the  first  weight.  The 
largest  coefficient  values  for  (B^J  are  in  the  range  of  +.16E14. 

When  the  weights  are  spread  out  over  a  wider  range,  larger 
values  of  M  can  be  tolerated  in  sum  (24),  without  encountering 
significant  round-off  error.  For  example,  a  set  of  N  »  16 
uniformly  distributed  random  weights,  over  the  (0,1)  interval, 
were  utilized  in  figure  15  without  any  problems.  But  when  M  was 
increased  to  20  in  figure  16,  again  for  uniformly  distributed 
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weights,  the  upper  edge  of  the  ROC,  for  >  .99,  was  useless. 
Nevertheless,  a  significant  portion  of  the  ROC  for  lower  ?D 
values  is  still  acceptable. 

The  lesson  to  be  drawn  from  these  results  is  that  the  partial 
fraction  expansion,  leading  to  the  exceedance  distribution 
function  in  (24),  has  utility  for  spread  out  weights  (wroJ  and 
moderately  low  values  of  M,  the  number  of  envelope-squared 

4 

samples.  However,  it  will  not  be  a  viable  tool  for  large  values 
of  M,  nor  for  general  weight  structures  which  may  have  some  close 
or  equal  values.  The  more  general  approach  presented  in  [2],  in 
£<STrms  of  an  arbitrary  characteristic  function,  has  no  such 
limitations,  on  the  other  hand,  although  the  numerical 
calculations  required  are  more  extensive. 
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Probability  of  False  fil arm 

Figure  12.  ROC  for  ,  r=.9S 
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Probability  of  False  fllarm 
Figure  13.  ROC  -for  M-8 ,  r=.39 
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Probability  of  False  Alarm 

Figure  14.  ROC  for  M=9 ,  r=.99 
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Probability  of  False  Alarm 
Figure  16.  ROC  -for  M-20,  Random  Weights 
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CHI-SQUARED  APPROXIMATION  FOR  ARBITRARY  WEIGHTS 

The  difficulty  of  evaluating  the  ROC  from  exact  character¬ 
istic  functions  of  the  form  of  (7)  and  (10)  has  prompted  the  use 
of  approximations  that  attempt  to  extract  an  effective  number  of 
independent  samples  from  a  general  weight  structure,  and  use  this 
parameter  in  a  simpler  chi-squared  fit.  For  example,  in 
(6;  (38)  and  sequel],  such  an  approximation  was  fruitfully 
employed  to  study  the  stability  of  a  spectral  analysis  technique 
employing  equi-weighted  overlapped  segments.  Also,  in 
{9;  (A-24)  -  (A-28)],  a  chi-squared  approximation  was  adopted  for 
the  analysis  of  a  diversity  combiner  in  a  partially-correlated 
fading  channel.  However,  in  this  latter  case,  no  quantitative 
measure  of  the  error  in  the  approximation  was  given. 

PARAMETERS  OF  APPROXIMATION 

Here,  we  will  address  the  adequacy  of  the  chi-squared 
approximation  for  a  general  exponential  weight  structure  of  the 
form  of  (27).  We  begin  by  generalizing  the  chi-squared 
characteristic  function  in  (12)  to  the  candidate  form 

-«e 

feU)  ■  (1  -  Uwea)  ,  (28) 

where  w^  is  an  effective  weight  and  M„  is  an  effective  number  of 
e  a  e 

envelope-squared  samples,  which  may  be  noninteger.  (The  number 
of  degrees  of  freedom  in  (28)  is  2Mg.)  The  corresponding 
probability  density  and  exceedance  distribution  functions  are 
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P«(u> 


N  -1 

u  exp 


(wea) 


M. 


for  u  >  0  , 


HMe)  (wea) 


Qe(u)  -  r[Me,  “)/r<Me)  for  u  >  0  , 


(29) 


respectively,  where  !*(•,*)  is  the  incomplete  gamma  function 
{10;  6.5.31.  These  results  generalize  (13)  and  (14).  The 
(scaled)  cumulants  of  this  gamma  distribution  follow  from  (28)  as 

TTrrrr  Vkl  -  Me  <"ea)k  £or  k  * 1  •  <30> 

The  mean  and  variance  of  this  approximation  are  therefore 

2  2 

H  w  a  and  M  w  a  ,  respectively. 

VC  V  V 

When  we  equate  these  first  two  moments  of  the  generalized 
chi-squared  approximation  <28)  to  the  first  two  moments  of 
decision  variable  x  in  (9)  and  (8),  we  find 


i 


w 


e 


M. 


3 

w. 


n 

2 _ 

.Dl*  1  y 

r>2 

L — ^  m 
m»  1 


(31) 


For  example,  if  all  the  weights  are  equal,  then  Mg  »  M.  On  the 
other  hand,  if  all  the  weights  are  zero  except  for  one,  then 
Mg  »  1 .  Both  of  these  limiting  cases  obviously  agree  with 
physical  intuition.  Observe  that  wg  and  are  independent  of 
parameter  a  or  R,  the  signal-to-noise  ratio. 

For  the  exponential  weight  structure  in  (27),  the  effective 
number  of  weights  and  the  effective  weight  are 
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M 


H. 


1  r  1  -  r 
1  -  r  i  ♦  ,» 


w  «  rj—  for  W. 
e  i 

e 


(32) 


It  should  be  noted  th«ic  as  H  ■»  »,  effective  number  saturates 
at  value  (1  +  r)/(l  -  r),  which  is  not  infinite. 

Since  the  incomplete  gamma  function  in  (29)  is  tedious  to 
compute  for  noninteger,  performance  could  be  bracketed  by  the 
two  cases  ,  K^+l,  where  is  the  integer  part  of  Me-  Or 
interpolation  could  be  used  between  these  two  cases.  Instead,  we 
shall  choose  examples  for  which  Mg  is  an  integer;  this  allows  us 
to  use  a  form  like  (14),  which  is  easily  computed  upon 
replacement  of  M  by  Mg. 


GRAPHICAL  RESULTS 

The  first  example  of  the  use  of  a  chi-squared  approximation, 
‘for  the  exponential  weight  structure  in  (27),  is  furnished  by 
figure  17  for  M  »  5,  r  ~  .69388907;  this  particular  r  value  is 
chosen  to  yield  M  -  4,  asmay  be  verified  from  (32).  The  exact 
results  (solid  lines)  in  this  figure  were  obtained  by  the  method 
of  the  previous  section,  namely,  all  weights  different.  The 
three  dashed  curves  arc  yielded  by  the  chi-squared  approximation 
of  this  section,  with  Mg  -  4;  the  latter  are  seen  to  be 
optimistic  by  almost  1  dB  along  the  left  edge  of  the  figure. 

When  M  is  increased  to  25  and  r  decreased  to  .60000182,  again 
resulting  in  M  •  4,  figure  18  shows  that  the  chi-squared 
approximation  is  far  worse.  The  reason  for  this  behavior  is  ^nat 
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25  significantly  different  weights  cannot  be  well  represented  by 
4  equal  weights  in  terms  of  evaluating  the  detection  capability 
of  the  energy  detector  (1). 

The  series  of  plots  in  figures  19,  20,  21,  22,  23  correspond, 
respectively,  to  Me  «  8,  16,  32,  64,  128,  for  various 
combinations  of  M  and  r,  as  indicated  on  the  figures.  Again,  the 
chi-squared  approximation  is  generally  optimistic  in  the  useful 
range  of  performance.  For  M  »  64  in  figure  20,  the  discrepancy 
is  almost  1  dB  along  the  left  edge.  However,  for  large  M,  like 
200  in  figure  23,  the  difference  is  only  about  .25  dB  along  the 
left  edge. 

The  results  in  figures  21,  22,  23  for  Me  -  32,  64,  128, 
respectively,  were  not  obtainable  from  the  all-weights-dif f erent 
method  of  the  previous  section,  due  to  excessively  1  irge 
coefficients  {B^}  in  (22).  Instead,  it  was  necessary  to  resort 
to  the  numerical  integration  procedure  given  in  [2J;  the  values 
of  increment  Aj-  and  length  Lg.  appropriate  to  each  case  are 
indicated  on  each  figure. 

A  conclusion  to  be  drawn  from  the  results  in  this  section  is 

that,  although  the  chi-squared  approximation  is  much  better  than 

the  Gaussian  approximation,  it  is  still  not  adequate  for  accurate 

performance  predictions  within  a  few  tenths  of  a  decibel.  The 

chi-squared  approximation  is  generally  unacceptable  for  small  Me, 

unless  r  is  very  close  to  1.  And  for  large  ,  it  is  acceptable 

in  some  regions  of  the  ROC,  but  not  in  others,  especially  if  the 

M-l 

extreme  weight  ratio,  r  ,  is  very  small. 
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Probability  of  False  Alarm 

Figure  17.  ROC  for  M*5 ,  r». 69388907  (Me-4) 
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Probability  of  Detection 
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Probability  of  False  fllarm 

Figure  19.  ROC  for  h=10,  r». 83623826  CMe=8) 
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Probability  of  False  Alarm 

Figure  20.  ROC  for  M=64,  r=. 88242683  (Me«lG) 
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Probability  of  False  Alarm 

Figure  21.  ROC  for  M=50f  r«. 94648071  (Me=32) 

39 


Probability  of  Detection 


TR  8753 


Probability  o-f  False  Alarm 
Figure  22.  ROC  for  M«100f  r*. 97288022  CMe=64) 
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Probability  of  False  Rlarm 
Figure  23.  ROC  for  M«200,  r-. 98634730  (Me-128) 
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THIRD-ORDER  APPROXIMATION  FOR  ARBITRARY  WEIGHTS 

When  a  constant  c  is  added  to  a  random  variable,  the 
characteristic  function  is  modified  by  multiplication  by  the 
factor  exp(ict)*  Accordingly,  a  further  generalization  of  the 
chi-squared  characteristic  function  in  (28)  is  afforded  by 

exp(i£b  a) 

fc(  l)  *  - jp  -  exp^i£bca  -  Mc  ln(l  -  i£wca)J  .  (33) 

(1  -  Hwca)  c 

This  form  now  has  three  parameters  to  choose,  namely  w  ,  b  ,  and 

V*  w 

effective  number  of  samples  Mc.  This  is  in  distinction  to  the 
chi-squared  approximation  (28)  and  the  Gaussian  approximation 
(A-2),  both  of  which  had  only  two  free  parameters  to  adjust. 

Thus,  whereas  we  only  matched  the  first  two  moments  in  (30)  and 
( A- 3 ) ,  respectively,  to  those  of  decision  variable  x,  we  can  now 
match  the  first  three  moments  of  x  if  we  use  characteristic 
function  model  (33). 

The  cumulants  of  characteristic  function  (33)  are 


xcU)  ■  Mc  wc  a  +  bc  a  ' 


nr^TJ  xc(k)  -  "c  (wc  a)k  £oc  k  2  2 


(34) 


When  the  first  three  cumulants  (or  moments)  of  (34)  are  equated 
with  the  corresponding  quantities  of  decision  variable  x,  as 
given  by  (8),  the  unique  solutions  for  the  parameters  of  (33)  are 


M 


c 


b 


c 


w 


1 


(35) 
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where 


(36) 


It  should  be  noted  that  the  parameters  in  (35)  are  independent  of 
parameter  a  or  R,  the  signal-to-noise  ratio. 

The  probability  density  function  corresponding  to 
characteristic  function  (33)  is 


"c-l 

(u  -  b  a)  exp 


u  +  b„aN 
c 


Pc(u) 


wca 


r("c)  (wc») 


M 


for  u  >  b  a  ,  (37) 

w 


and  zero  otherwise.  The  exceedance  (gamma)  distribution  function 
is  an  obvious  generalization  of  (29) ,  or  (14)  if  Mc  is  integer; 
see  [10;  6.5.3,.  6.5.2,  6.5.13]. 


Qc(u) 


V 


u  -  bca>i 


wca 


n«c> 


(u  -  b^' 


'«c-l 


wca 


for  u  >  b  a  .  { 38 ) 
c 


For  threshold  value  T,  the  false  alarm  and  detection 
probabilities  follow  immediately  as 


'"c-l 


(T  -  b^'l 

fT  -  b  a^ 

c 

p  •  p 

c 

w 

'  *D  M  -1 

w  a 

k  C  / 

c 

l  c 

(  39) 


provided  that  T  >  bca. 
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EXPONENTIAL  WEIGHTS 

We  now  restrict  attention  to  the  exponential  weight  structure 
w  ■  l — ~  for  1  <  m  <  M  ,  with  t  *  rM  ,  (40)  . 

m  i  *  l  ~ 


where  we  have  normalized  at  W^  »  1.  Then,  from  (36)., 


W, 


U  t  t  t2  + 
1  +  r  +  r2  + 


+  t 


k-1 


+  r 


k-1 


(41) 


In  particular. 


W, 


W-.  - 


1  -  r  1  +  t 

1  -  t  1  +  r 


w. 


a  -  n2  i  +  t  +  t2 
“  I r I  - J  *  (42) 

J  1  +  r  + 


The  parameters  in  (35)  then  follow  by  substitution  as 


N 


(l  -  r3)2  (l  -  t2)3 


(l  -  r2)  u 


-  t2)3  .  1  -  t  a  ♦  f|3  2 

-  t3)2  1  '  ‘  “  +  '  U  +  t  +  t2J 


(43) 


w 


1  -  r2  1  +  t  +  t2 

- - T  2  ' 

1  -  t*  1  +  r  +  t4 


b_  « 


(1  -  r  t)  (r  -  t) 

(1  ♦  r)2  (l  +  t  +  t2) 


(44) 


For  equal  weights,  wffl  -  1/M,  we  get  the  usual  reduction  to 
W^  ■  1,  Wj  •  1/M,  W^  -  1/M2,  giving  Mc  »  M,  w  *  1/M,  bc  *  0. 
Furthermore,  it  is  shown  in  appendix  B  that  additive  constant  bc 
in  (33)  and  (37),  as  determined  from  (35)  and  (36),  is  never 
negative,  for  any  nonnegative  weight  structure  . 
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GRAPHICAL  RESULTS 

The  first  example  we  consider  here  is  M  *  25,  r  -  .75049209, 
for  which  (43)  gives  M  ■»  4 ;  again,  the  reason  for  the  particular 
choice  of  r  is  made  so  that  Mc  is  integer  and  (39)  can  be  used. 
The  approximation  afforded  by  (39)  is  superposed  (dashed  lines) 
in  figure  24  on  the  exact  results  (solid  lines)  obtained  from 
(25).  Increasing  M  to  64  and  changing  r  to  .75049170,  so  that  M 

w 

is  maintained  at  4,  generates  virtually  the  same  approximation. 

The  fit  is  poor  and  rather  optimistic  at  the  left  edge  of  the 

figure,  due  to  the  small  value  of  Mc,  namely  4. 

For  M  «  50  and  r  «  .96915298,  M  is  increased  to  32  and  the 

results  are  compared  in  figure  25.  Now,  the  fit  afforded  by  the 

constant  plus  chi-squared  approximation  is  rather  good  over  the 

entire  range  of  false  alarm  and  detection  probabilities  shown,*  in 

fact,  the  approximation  is  optimistic  by  about  .1  dB  on  the  left 

edge  of  the  figure.  The  reason  for  this  development  is  the 

larger  value  of  the  effective  number  of  samples,  Mc,  namely  32. 

Two  more  results,  for  Mc  equal  to  64  and  128,  yield  similar 

conclusions  in  figures  26  and  27,  respectively.  Again,  the 

exponential  weight  structure  was  employed.  However,  the  goodness 

of  fit  of  the  constant  plus  chi-squared  approximation  is  not 

limited  to  this  type  of  weights,  but  in  fact  applies  to  arbitrary 

structures.  To  back  up  this  statement,  an  example  of  uniformly 

distributed  random  weights  for  M  ■  133  and  M  -  77.971  is 

c 

displayed  in  figure  28;  the  overlay,  which  used  Mc  ■  78  in 
approximation  (39),  is  seen  to  be  very  good  for  this  value  of  Mc. 
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Probability  of  False  fll arm 
Figure  24..  ROC  for  M=25,  r*. 75049209  (Mc=4) 
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Probability  of  False  Alarm 
Figure  25.  ROC  for  M=50,  r=. 96915298  (Nc=32) 
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Probability  of  Detection 


Probability  of  False  Rlarm 

Figure  26.  ROC  for  M*100,  r=. 98445939  (Mc=64) 
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Probability  of  False  Alarm 

Figure  27.  ROC  for  M«200,  r=. 99220012  CMc*l28) 
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Figure  28.  ROC  for  M=133,  Random  Weights  (Mc=?8) 
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APPLICATION  TO  EIGENVALUE  PROBLEH 


Earlier,  in  (10)  and  (8;  (24)],  a  particular  characteristic 
function  was  given  which  has  occurred  in  a  number  of  statistical 
analyses.  That  characteristic  function,  in  normalized  form,  is 


£XU> 


(45) 


where  R  is  the  per-sample  signal-to-noise  ratio  and  {X^}  are  the 
eigenvalues  of  the  normalized  covariance  matrix  P  of  the  fading 
signal  components.  By  expanding  the  In  of  (45)  in  a  power 
series  in  it,  the  cumulants  of  random  variable  x  are  found  to  be 


M 

( k-1 ) 1  Xx<k)  ” 


M  k 


(1  +  R  X  « 
m 


m»l 


HE  o *n^- 


m»l  n«*0 


*  h  +  s"  tr(pn)  for  k  i 1  - 


(46) 


n«l 


where  we  have  used  the  simplifying  result  in  appendix  C  regarding 
sums  of  powers  of  eigenvalues,  in  particular,  there  follows 
from  (46),  the  first  three  cumulants  of  x  in  terms  of  tr(Pn): 


xx(i> 

-  M 

+ 

R 

t  r  ( P ) 

9 

Xx(2) 

-  M 

+ 

2R 

tr(P) 

*  R2 

tr ( P2  )  , 

k<3) 

-  M 

+ 

3R 

tr(P) 

+  3R2 

tr  ( P2  )  +  R3  t  r  (  P 3  )  . 

(47) 
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PARAMETERS  FOR  CANDIDATE  APPROXIMATION 

In  this  section,  we  will  approximate  exact  characteristic 
function  (45)  by  the  form  employing  the  constant  plus  chi-squared 
idea  again,  namely 

exp{ i£bj)  ,  \ 

fdU)  ,  - _  exp[i^bd  -  Md  ln(l  -  Uwd)J  .  (48) 

(1  -  Hwd)  d 

The  cumulants  are  given  by  a  form  very  similar  to  (34),  and  in 
particular,  the  first  three  (scaled)  cumulants  of  characteristic 
function  (48)  are 

Xd(l)  “  Md  wd  +  bd  '  xd*2*  "  Md  wd  '  §xd(3)  “  Md  wd  ’ 

If  the  first  three  cumulants,  XdUO  for  k«l,2,3,  were  specified, 
we  could  then  solve  (49)  for  the  required  parameters  according  to 


Now,  we  set  the  cumulants  of  approximation  (48)  equal  to  the 
exact  cumulants  given  by  (47),  and  then  solve  (50)  for  the 
required  parameter  values.  Then,  approximation  (48)  to  exact 
characteristic  function  (45)  is  available  for  numerical 
evaluation.  If  cumulants  (xx(X)}  for  k-1,2,3  can  be  evaluated 
either  analytically  (via  eigenvalues  in  (46)  or  by  the  trace 

relations  in  <47))  or  numerically  (estimated  via  finite  time 
averages),  then  the  parameters  in  (50)  can  be  determined  and  the 
corresponding  ROC  found. 


*du> 

(  3  )/2  j 


Xd(3)/2 

xdU) 


xd(D 


xd<2> 

Xd<3)/2 


(50) 
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EXACT  PERFORMANCE  OF  (45) 


If  signal-to-noise  ratio  R  *  0  in  (45),  then 

—  M 

f  (fc)  »  (X  -  i£)  and  there  follows,  in  a  manner  similar  to 
(14),  Qx(u)  *  Em_i(u)  and  Pp  »  EM_1(T)  for  threshold  T. 

If  R  >  0  and  all  the  eigenvalues  {X^}  in  (45)  are  distinct, 
then,  in  a  manner  similar  to  (21),  we  can  express 


fx(U 


M 


y  VR) 

fcl  1  -  + 


t(l  +  F  Xj  ’ 


(51) 


where  coefficients 


VR) 


l1  *  R  Xm) 


M-l 


M 

T1  (Xm  -  Xk) 


for  1  <  m  <  M 


(52) 


k-1 

k*m 


The  exceedance  distribution  function  is  then 


Qx(u) 


M 

Evr)  e>P (l  T  B  X  ) 
m-l 


for  u  >  0 


R  >  0  , 


(53) 


and  the  detection  probability  is 

M 

pd  ■  Ev*>  e,p(r Hi)  foc  T  >  0  '  R  >  0  •  l54> 

m-l 

The  eigenvalues  {X^}  of  normalized  covariance  matrix  P  are 
independent  of  signai-to-noise  ratio  R;  however,  coefficients 
(Bm(R))  are  dependent  on  R  and  explicitly  indicated  so. 
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GRAPHICAL  RESULTS 

The  only  example  that  we  consider  here  is  a  covariance  matrix 

P  -  (p _ 1,  where  pm_  «  plm~nl.  in  particular,  for  M  -  10  and 

sin  in 

p  >  .5,  the  M  eigenvalues  {Xffl}  of  P  were  evaluated  and  the 
results  on  page  55  were  used  for  an  exact  evaluation  of  the 
detection  and  false  alarm  probabilities;  these  are  displayed  as 
solid  lines  in  figure  29. 

Then,  we  returned  to  matrix  P,  ignored  the  knowledge  of  the 
eigenvalues,  and  instead  employed  the  trace  relations  in  (47)  and 
appendix  C  to  evaluate  the  cumulants  of  random  variable  x.  These 
were  substituted  in  (50)  to  determine  the  parameters  of 
characteristic  function  (48),  as  explained  in  the  sequel  to  (50). 
Then,  the  method  of  (2)  was  used  to  obtain  the  corresponding  ROC. 

These  results  are  overlaid  as  dashed  lines  in  figure  29,  for 
three  selected  values  of  signal-to-noise  ratio  R  (in  decibels). 
The  agreement  for  small  signal-to-noise  ratios  is  very  good,  and 
can  be  explained  by  observing  that  (45)  approaches  the  chi- 
squared  characteristic  function  in  this  case.  Approximation  (48) 
is  also  excellent  for  very  small  false  alarm  probabilities, 
despite  the  fact  that  the  equivalent  number  of  samples,  Md,  is 
rather  small;  for  example,  the  three  curves  in  figure  29  for 
R  «  2,5,8  dB  have  Md  -  5.79,  4.83,  4.31,  respectively. 

Another  example  for  M  ■  32,  p  *  .5  is  displayed  in  figure  30. 
Here,  the  values  of  Md  for  the  four  overlays,  R  -  -2, 0,2, 4  dB  are 
24.1,  20.5,  17.6,  15.4,  respectively.  These  larger  values  of  Md 
account  for  the  improved  fit  to  the  exact  results. 
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FOURTH-ORDER  APPROXIMATIONS  FOR  ARBITRARY  WEIGHTS 

In  this  section,  we  will  consider  a  couple  of  fourth-order 
fits  to  a  specified  characteristic  function  and  will  match 
cumulants  (or  moments)  through  fourth-order. 


GAUSSIAN  PLUS  CHI-SQUARED  FIT 


The  initial  fourth-order  fit  of  interest  here  corresponds  to 
the  characteristic  function  of  a  (nonzero  mean)  Gaussian  random 
variable  plus  a  chi-squared  variate.  That  is,  the  candidate  is 


ffm 


exp  fi  (b,  —  xt  c  .1  /  15  n 

- --  -f— -  e*p(ub£  -  £t2cf  -  Mf  InU  -  itw£))  . 

(l  -  Uw£) 


(55) 


The  first  four  cumulants  of  characteristic  function  (55)  are 


Xf(l)  -  bf  +  Mf  wf  ,  XfU)  -  cf  +  Mf  wf  , 


|xf(3)  -  Mf  wjj  ,  |xf(4)  -  Mf  wj  . 


(56) 


If  the  cumulants  are  specified,  the  parameters  for 
characteristic  function  (55)  can  be  determined  explicitly  as 


Mf  - 


bf  -  X«(D  - 


(xf(3)/2) 

(xf(4)/6) 

Jx£<3)/2V 


Xf(4)/6 

I  '  wf  "  xTTTTTS  ' 


f'*'  '  (x£H>/6)2  '  “  X‘(2’  _  X£M>/6 


(xf(3)/2j 


(57) 


Numerical  results  will  be  presented  in  a  later  section, 
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NON-CENTRAL  CHI-SQUARED  FIT 


The  other  fourth-order  fit  that  we  consider  corresponds  to  a 
generalized  non-central  chi-squared  variate,  namely 
characteristic  function 


Ub, 


exp 


1  -  ite. 


,  \M« 
(l  -  1ft.,)  9 


exp 


Ub, 


-  "9  lnU  -  UV 


(58) 


This  is  called  generalized  because  we  do  not  force  c^  •  w^. 
The  In  of  (58)  can  be  expanded  in  a  power  series  in  it: 


+«>  +«> 

in  yt)  -  «bg  uyl  *  Ng  ^|utvg)k  . 

j-0  k-1 


(59) 


The  first  four  cumulants  of  this  characteristic  function  are  then 


xg(l)  -  bg  .  «g  »g  ,  Xg(  2 )  -  2  bg  cg  ♦  »g  w‘  , 


ix„(3)  -  3  b  c*  +  K  w*  lx„H)  -  4  b„  c9  +  w*  (60) 

2g  gg  gg  6g  gg  gg 


The  inversion  of  these  nonlinear  equations,  for  the  parameters  in 
terms  of  the  cumulants,  is  not  possible  in  closed  form,  as  it  was 
for  candidate  characteristic  function  (55).  This  limitation 
tends  to  discourage  use  of  the  non-central  chi-squared 
approximation  (58).  However,  in  appendix  D,  an  efficient 
numerical  procedure  for  solving  (60)  for  the  required  parameters 
is  developed  and  programmed.  Application  of  this  approximation 
procedure  is  deferred  to  a  later  section. 
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PERFORMANCE  IN  STEADY  STATE  NOISE 

Up  to  this  point,  the  number  of  samples,  H,  has  been  finite, 
both  for  signal-present  as  well  as  signal-absent?  then,  the  noise 
output  of  the  exponential  integrator,  (27}  or  (40),  has  not 
reached  steady-state.  In  this  section,  the  number  M  of  noise 
samples  will  be  set  equal  to  «,  thereby  allowing  the  integrator 
noise  output  to  reach  steady  state.  However,  the  number,  N,  of 
samples  containing  signal  (if  present}  will  remain  finite. 

This  situation  arises  in  practice,  for  example,  when  the 
precise  arrival  time  of  the  signal  is  unknown.  The  use  of 
surplus  envelope-squared  samples  for  m  >  N,  does  not 

improve  performance,  since  these  particular  samples  are  always 
noise-only;  in  fact,  these  extra  samples  always  degrade 
performance,  the  exact  amount  depending  on  the  relative  sizes  of 
weights  (wm)  for  m  >  N  compared  to  m  <  N.  Here,  we  will  give  a 
method  for  quantitatively  assessing  the  impact  of  these  surplus 
noise-only  samples  on  the  operating  characteristics. 
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CHARACTERISTIC  FUNCTION 


The  characteristic  function  of  the  decision  variable  is  an 
obvious  generalization  of  (7)  to  the  form 


- 


ft l1  -  i(v.) 


-1 


(61) 


where  the  signal-to-noise  ratio  parameter  aB  now  takes  the  form 


a 


m 


{  1 
U  +  R 


for  noise-alone 


r  for  signal-plus-noise> 

IB 


for  1  <  m  <  M  « ■  *  .  (62) 


The  particular  case  that  will  be  considered  at  length,  here,  is 
that  of  a  finite-duration  constant-strength  signal,  which  is 
accommodated  mathematically  by  setting 


fR  for  1  <  ra  <  N  'j 
m  10  for  N  <  m  <  M  - 


(63) 


When  signal-to-noise  ratio  R  is  equal  to  zero,  that  is, 
signal-absent,  the  characteristic  function  in  (61)  reduces  to 


?XU) 


,-l 


n<l '  Uw«) 


lra-1 


(64) 


Unfortunately,  even  for  the  exponential  averager. 


wm  -  (1-r)  r 


m-1 


for  1  <  m  <  M  -  •  , 


(65) 


the  noise-only  characteristic  function  in  (64)  takes  a  form. 


62 


7R  8753 


lxU)  « 

which  is  not  expressible  in  closed  focm;  see  [11;  (89.18.3)]. 
(Likewise,  the  finite  product  cannot  be  simplified;  see 
[11;  (89.18.2)].)  This  necessitates  termination  of  the  infinite 
product  in  (66),  being  sure  'to  keep  the  remainder  below  an 
acceptable  tolerance;  this  issue  iu  addressed  in  appendix  E. 


.m*i 


T  **■ 


w 

no- 


i K  ( 1-r  )  r 


m-1' 


(66) 


CUMULANTS 


For  general  characteristic  function  (61),  the  cumulants  are 


(k-l)l 


X„(k) 


E 

m-l 


(w  *0  for  1  <  k 
mm  ** 


(67) 


For  the  special  case  of  the  exponential  averager  (65)  and  the 
finite-duration  signal  (63),  these  cumulants  reduce  to 


(k-l)l  Xx(k) 


(1  -  r) 


1-r 


(1  +  R)k(l  -  rkN]  +  r 


kN 


(68) 


At  the  same  time,  characteristic  function  (61)  becomes 


£x:t» 


N 


no 

m-l 


iE.  (1-r)  rm-1(l  +  R) 


w 

ini- 


if.  (1-r)  r 


m-l' 


m-N+1 


-1 


(69) 
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In  particular,  for  noise-alone,  then  R  -  0  and  (58)  reduces 
to 

_ i _  z  (j<)  .  Lk z Li _  - LA. i — £ ) -  (70) 

(k-l)»  k  k-1  * 

l)'  x  1  -  rK  1  +  t  +  ...  +  rK  1 


The  three  lowest-order  cases  are 


Xv(l)  -  1  ,  Xv<2)  -  '  yXJ3)  " 

x  x  i  +  r  &  x  l+r  +  rz 


(71) 


For  signal-present,  R  >  0,  the  three  lowest  cumulants  are, 
from  (68), 


Xx(l)  -  1  +  R  -  R  rN  , 


V2)  "  rT~f[{1  +  r)2(1  ■  r2N)  +  r2N]  * 


ixx(3)  -  J-1  +  R)3fi  _  r3N]  +  r3N]  .  (72) 

,4X  1  +  r  +  r/L  k  '  J 


Here,  N  is  the  number  of  signal  components,  R  is  the  signal-to- 
noise  ratio  per  sample,  and  r  is  the  exponential  decay  factor  for 
the  weight  structure  (65). 

In  the  evaluation  of  the  signal-present  characteristic 
function  (69),  the  second  product  will  have  to  be  terminated  at  a 
finite  limit  m  »  L  (>  N).  The  error  due  to  this  truncation  is 
addressed  in  appendix  E. 
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GRAPHICAL  RESULTS 

An  example  of  the  results  in  this  section  for  M  «  N  »  32, 
r  »  .9,  is  displayed  in  figure  31,  as  obtained  via  exact  results 
(66)  and  (69),  along  with  the  truncation  procedure  of  appendix  E. 
Superposed  as  dashed  lines  are  the  results  of  using  the  constant 
plus  chi-squared  approximation  (48),  where  the  parameters  are 
ol  ained  from  the  cumulants,  according  to  (50).  The  cumulants 
themselves  are  given  by  (72).  The  effective  number  of  samples, 
in  (48),  takes  on  the  values  10.680,  10.676,  10.673,  and 
10.672  for  the  four  signal-to-noise  ratios  of  0,  2,  4,  and  6  dB 
indicated  in  the  figure.  This  relatively  small  value  of  is 
the  reason  for  the  discrepancy  in  figure  31  between  the  exact  and 
approximate  results. 

Figure  32  is  drawn  for  M  »  •,  N  »  50,  and  r  -  .96915298; 
compare  figure  25,  for  which  Mc  »  32.  The  values  of  Md  for  the 
three  overlaid  curves  for  signal-to-noise  ratios  equal  to  -2,  0, 
and  2  dB  are  33.531,  33.030,  and  32.624,  respectively.  These 
larger  values,  for  the  effective  number  of  samples,  lead  to 
better  agreement  in  this  figure;  in  fact,  the  approximation  is  in 
error  by  only  .15  dB  along  the  left  edge  of  the  figure. 
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Probability  of  False  Alarm 
Figure  32.  ROC  for  M=<»  N=50,  r=. 96915298 
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BLOCK  EXPONENTIAL  WEIGHTING 

% 

In  this  section,  we  again  consider  a  weighted  energy  detector 
in  steady  state,  that  is,  M  *  «*.  However,  the  averager  now 
operates  on  blocks  of  data  points  which  are  equally  weighted,  but 
which  are  themselves  exponentially  weighted.  That  is,  the 
decision  variable  x  is  now  given  by 


where  the  weights  {wm)  are 


for 

for 

for 


1  <  m  <  b 
B  <  m  <  2B 
2B  <  m  <  3B 


(73) 


(74) 


Here,  B  is  the  block  size  and  the  weights  have  been  normalized 
at  Wj  «  1.  The  following  diagram  illustrates  the  block 
exponential  weighting  structure. 
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SIGNAL  STATISTICS 

The  signal,  if  present,  occupies  the  first  N  samples  of  sum 
(73),  where 

3  ■  |  (75) 

is  presumed  integer;  that  is,  3  is  the  number  of  blocks  occupied 
by  signal  (when  present).  The  signal-to-noise  ratio  parameter  is 

'  1  for  noise-only  " 

a^  »  <  for  1  <  m  <  N  ,  (76A) 

J.+R  for  signal  plus  noise. 

and 

a^  «  1  for  N  <  m  <  *  .  (76B) 

m 

CHARACTERISTIC  FUNCTION 

The  characteristic  function  of  x  in  (73)  for  signal  present 
is,  using  the  independence  of  the 

T1 

1  -  Uw»am)  - 

fit1  -  *«  ¥  r3  n  ♦  R>)  f[(l  -  lt  ¥  '3) 

j-0  3-J 

Again,  an  infinite  product  is  required  and  the  truncation 
procedure  given  in  appendix  E  is  directly  relevant. 
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CUMULANTS 


The  cumulants  of  decision  variable  x  follow  readily  from 
(11),  upon  expansion  of  In  fx(£)  in  a  power  series  in  i£: 


,k-l 


( k-1 )  I 


Xx(k) 


LLlrXl/B)  ,  [(1+R)k  fl-rkJ1  +  rkJ]  for  k  >  1.  (78) 
l+r+- - L  v  '  J 


+  r 


The  four  lowest-order  cumulants  will  be  used  in  fitting  the  exact 
characteristic  function  (77)  by  approximations  (55)  and  (58). 


GRAPHICAL  RESULTS 

Results  for  the  operating  characteristics  for  B  »  4,  J  *  32, 
and  r  «  .95  are  presented  in  figure  33.  Thus,  from  (75),  the 
signal  (when  present)  occurs  on  N  •  128  samples.  Superposed  as 
dashed  lines  is  the  approximation  afforded  by  third-order  fit 
(33)  and  (39).  The  discrepancy  is  only  .1  dB  along  the  left  edge 
of  the  figure. 

Another  example  of  block  exponential  weighting,  for  B  ■  4, 

J  »  16,  and  r  «  .9,  is  displayed  in  figure  34.  The  dashed 
overlay  is  again  the  third-order  approximation  (33),  which  is 
optimistic  by  about  .15  dB  along  the  left  edge  of  the  figure. 

The  exact  results  from  figure  34  are  repeated  in  figure  35, 
but  now  the  overlays  are  the  two  fourth-cr'  r  approximations  (55) 
and  (58).  The  latter  two  approximat ionr,  are  indistinguishable 
from  each  other  over  the  entire  range  of  c  babilities  displayed. 
Furthermore,  they  differ  from  the  exact  results  only  by  .05  dB  at 
the  left  edge  of  the  figure. 
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Probability  of  False  Rlarm 
Figure  33.  ROC  for  B=4 ,  J=32,  r=.95 
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Probability  of  False  Rlarm 

Figure  35,  ROC  for  8=4,  J“16,  r=.9,  f ourtn-or der 
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SUMMARY 

The  receiver  operating  characteristics  of  a  variety  of 
weighted  energy  detectors,  for  Gaussian  signals  in  noise,  have 
been  investigated  exactly  and  compared  with  five  different 
approximate  procedures.  The  Gaussian  and  chi-squared 
approximations  have  been  found  to  be  generally  inadequate  for 
very  small  false  alarm  probabilities,  while  the  generalized 
chi-squared  (gamma)  and  both  fourth-order  fits  have  yielded  very 
good  results  over  the  entire  range  of  detection  and  false  alarm 
probabilities  considered.  The  only  limitation  of  the  latter 
approaches  is  the  need  to  have  additional  cumulants  (or  moments), 
since  the  first  two  cumulants  are  not  always  entirely  adequate 
for  accurate  performance  predictions. 

If  the  exact  characteristic  function  for  the  decision 
variable  of  a  system  can  be  determined,  either  analytically  or 
numerically,  then  the  receiver  operating  characteristics  can  be 
accurately  evaluated  by  the  method  of  [2J,  as  done  here. 

However,  there  are  occasions  where  it  may  be  desirable  or 
imperative  to  use  an  approximate  characteristic  function,  as  for 
example,  when  only  a  few  low-order  moments  are  known.  In  this 
fashion,  we  can,  for  example,  avoid  the  determination  of 
eigenvalues  or  avoid  the  evaluation  of  infinite  products.  Also, 
the  approximate  forms  will  frequently  be  faster  to  compute  than 
the  exact  results.  This  report  indicates  the  relative  accuracies 
inherent  in  some  of  the  standard  approximations  and  some  of  their 
generalizations,  which  should  be  considered  for  future  use. 
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APPENDIX  A  -  GAUSSIAN  APPROXIMATION 


The  characteristic  function  of  interest  was  presented  in 

M  H  ',_1 

fx(U  -  E{exp(  i  £x) }  -  n  w»  -  nc1  -  uw«a) 

m-1  .»*! 


where  {wmJ ,  for  1  <  a  <  M,  are  an  arbitrary  set  of  weights, 
aean  and  variance  of  random  variable  x  were  given  in  (9). 

Now,  if  energy  detector  output  x  in  (1)  were  a  Gaussian 
random  variable,  its  probability  density  function  would  be 


Pg<u> 


1 

(2n)  V 


exp 


/  2^ 
<*»  -  *q> 


2 


for  all  u  , 


where,  from  (9)  and  (4),  we  set 


1 

f  w  \ 

2 

2 

or 

or 

k  1  +  R  . 

.  (1  +  R)2  W2  , 

>  . 


The  exceedance  distribution  function  corresponding  to  (A-2) 


1* 

ffj  -  u') 

Qg ( u )  «  J  dt  Pg(t)  -  ♦ 

tl 

_ 

k  - 

for  all  u  , 


where 


${ t )  s  |  dv  exp(-v2/2) 

—  • 

is  the  normalized  Gaussian  cumulative  distribution  function 


(7): 

{ A-l ) 

The 

(A-2) 

(A-3) 

is 

( A-4 ) 

(A- 5) 
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At  this  point,  it  is  convenient  to  define  an  effective  number 
of  samples,  M0,  for  an  arbitrary  set  of  weights  {w^}  as  in  (31) 


Here,  we  used  (8)  and  (2). 

If  threshold  T  is  utilised  for  a  comparison  with  energy 
detector  output  x  for  a  decision  on  signal  presence  or  absence, 
then  the  approximate  false  alarm  probability  follows  from  (A-4): 

Pr  -  Qg<Tf  R-0)  -  #(m^<1-T))  ,  ( A- 7 ) 

with  the  help  of  (A-3)  and  (A-6).  Similarly,  the  approximate 
detection  probability  is 

?D  -  Qg(T;  R*0)  -  .(«J(l  -  tfj))  •  (A-8) 

Equations  (A-7)  and  (A-8)  produce  the  Gaussian  approximation  to 
the  operating  characteristics  of  the  energy  detector  (1), 
described  by  characteristic  function  (A-l).  They  depend  only  on 
the  single  parameter  Me  defined  in  (A-6),  in  addition  to  the  per- 
sample  signal-to-noise  ratio  R.  That  is,  M  and  (wm)  are  all 
collapsed  into  the  single  parameter,  effective  number  Me. 

An  immediate  obvious  problem  with  (A-8)  is  that  the  limit  of 
detection  probability  P&,  as  R  •*  «*,  is  not  1;  in  fact,  it  is 

<  1.  This  drawback  serves  as  a  warning  about  the  adequacy 
of  the  Gaussian  approximation. 
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For  the  approximations  in  (A-7)  and  (A-8),  we  can  explicitly 
solve  for  in  terms  of  Pp,  as  follows.  Let  *  be  the  inverse 
function  to  ♦;  see  [10;  26.2.23].  Then  (A-7)  can  be  solved  for 
threshold  T  according  to 

T  «  1  -  M *(PF)  .  ( A-9 ) 

Substitution  of  this  result  into  (A-8)  yields 


"e  R  *  !(pr)' 


1  +  R 


( A-10 ) 


It  now  follows  immediately  from  (A-10)  that,  for  specified 
and  PD,  the  required  signal-to-noise  ratio  R  is 


where 


(A— 11 ) 


F  *  ♦(Pp)  ,  D  ■  *<PD)  .  (A-12) 

The  result  in  (A-ll)  is  a  generalization  of  (1;  ( C-8 )  and  (11)] 
to  the  case  of  arbitrary  weights  {w^}.  It  is  immediately  obvious 
from  the  denominator  of  (A-ll)  that  the  desired  PD  must  be 
smaller  than 
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APPENDIX  B  -  POSITIVITY  OF  PARAMETER  fc>c 

Here,  we  will  show  that  the  parameter  bc  in  (35)  is  never 

negative,  regardless  of  the  weight  structure  {wmJ,  provided  that 

w  >0.  The  Cauchy-Schwartz  inequality  states  that 
m  ~ 


M  ^2  M  M 

E a- v  j  E  a»  E b” 

km»l  m-1  m-1 


( B— 1 ) 


for  any  real  quantities  {a^}  and  {bm).  If  we  let  am  «  wm 


3/2 


and 


b  -  »  1/,2(  then  (B-l)  yields 


M  2  M  M 

Un-1  m-1  m-1 


(B-2) 


that  is,  Wj  <  w3  w1#  where 


M 


Wt 


z_ ,  w»  • 

m-1 


(B-3) 


Therefore 


W2 

W1  “  4  ° 


( B-4 ) 


in  addition,  there  follows 


M  w 

c  c 


w2 

tq  <  w! 


( B-5 ) 
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APPENDIX  C  -  TRACE  RELATIONS  FOR  EIGENVALUES 

Suppose  MxM  matrix  P  -  [pmnl  has  eigenvalues  {\m},  1  <  m  <  M. 
Let  A  be  the  diagonal  matrix  of  eigenvalues  {A^}  and  let  Q  be  the 
normalized  modal  matrix  of  eigenvectors  of  P;  see  [12;  section 
1.13J .  Then  we  can  express  matrix  P  in  the  form 

P  -  Q  A  QT  ,  ( C-l ) 

from  which  there  follows  the  k-th  power 

Pk  «  Q  Ak  QT  .  ( C-2 ) 

He  now  use  the  trace  relation 


tr (A  B  C)  -  tr (B  C  A)  , 


( C-3 ) 


to  evaluate  the  trace  of  P: 


tr  (pk)  -  tr  (q  Ak  qt)  «  tr  (aR  QT  q)  -  tr(Ak)  "  •  (C-4) 


m«l 


That  is,  the  sum  of  the  k-th  powers  of  eigenvalues  { Xm }  can  be 

k 

obtained  from  the  trace  of  matrix  P  ,  without  ever  having  to 
evaluate  the  eigenvalues  at  all.  In  particular. 


M 


N 


-  tr(P)  -  ^„mn  . 
m-1  m«l 


( C-5 ) 


M  H 

- tc(p2)  -  £ » 


mn  nm 


( C-6 ) 


m- 1 


m,n-l 
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M  M 

E*»  ■  tt(p3)  ‘  E  P«n  Pnk  pkm  • 

m-1  m,n,k«l 

In  order  to  compute  the  sums  o£  the  three  lowest  powers  of  the 

eigenvalues  of  matrix  P,  we  simply  have  to  compute  the  three  sums 

on  the  elements  of  matrix  P,  as  indicated  in  (05)  through  (07). 

In  fact#  there  is  no  need  to  compute  matrices  P  or  P  either. 

* 

Thus#  a  seemingly  difficult  numerical  chore  is  replaced  by 
straightforward  simple  summations  of  products  of  matrix  elements, 
yielding  a  very  significant  savings  in  complexity  and  time. 
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APPENDIX  D  -  INVERSION  OF  EQUATION  (60) 


For  notational  efficiency,  we  suppress  all  the  g  subscripts 
in  (60),  let  «  x(k)/(k-l)S,  and  set  p  »  M  w.  The  nonlinear 
equations  then  take  the  form 


n  "  b  +  P  <  y2-2bc+pw, 
y,  «  3  b  c2  +  p  w2  ,  y,  ■  4  b  c3  +  p  w3  .  {D-D 

We  solve  the  first  two  equations  for  p  and  b,  getting 


y2  -  2  ?i  c 

w  -  2  c 


b 


*i  *  :  y2 

w  -  2  c  '  * 


(D-2) 


These  quantities  are  nov  substituted  in  the  third  and  fourth 
equations  in  (D-l),  resulting  in  the  highly  nonlinear  pair  of 
coupled  equations  for  c  and  w: 


c2  3  (yx  w  -  y2)  +  c  2  (y3  -  yx  w2)  w  <y2  w  -  y^)  -  0  ,  (D-3) 

3  3  2 

c  4  (y1  w  -  y2)  +  c  2  (y4  -  yx  w  )  •  w  (y2  w  -  y4)  -  0  .  (D-4) 


The  procedure  we  have  adopted  for  solving  these  latter  two 
equations  is  to  start  with  an  initial  guess  for  w  as  in  {57), 
namely 


w 


XUlZg 

X'3)/2 


(D-5) 


then  solve  quadratic  (D-3)  for  c;  substitute  this  result  into 
(D-4)  and  compute  the  left-hand  side;  now  vary  w  until  the 
left-hand  side  equals  zero.  Repeat  these  operations  until  c  and  w 
stabilize.  Equation  (D-2)  can  now  be  used  to  get  final  values  of 
p  and  b.  This  ib  the  numerical  procedure  used  in  the  main  text. 
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APPENDIX  E  -  TERMINATION  OF  INFINITE  PRODUCT 

If  we  terminate  the  infinite  product  for  the  characteristic 
functions  in  (66)  or  (69)  at  limit  value  m  •  Z.  (>  N),  then  the 
neglected  remainder  product  in  the  denominator  is 

Rem  -  1~Jfl  ~  i£  (1-r)  rm“1)  «  1  -  i£rL  -  £2  +  o(r3L].(E-l) 

m»L+l 


This  relation  enables  a  choice  of  L  to  control  the  neglected 

remainder.  For  example,  £  «  200,  L  «  220,  r  «  .9  leads  to 

2 

Rem  «  1  -  il.7E-8  -  1.4E-16.  Thus,  the  £  term  and  above  can  be 
safely  ignored.  One  final  product  in  the  denominator  of  (66),  by 
the  factor  1  -  i£rL,  will  account  for  Rem  and  suffice  for 
complete  accuracy,  up  to  computer  round-off  error  in  the 
characteristic  function  evaluation.  For  larger  values  of  r,  it 
is  necessary  to  increase  the  limit  L;  for  example,  £  -  150, 

L  -  700,  r  -  .96915298  yields  Rem  -  1  -  i4.5E-8  -  IE-15. 

If  we  terminate  the  infinite  product  for  the  characteristic 
function  in  (77)  at  limit  value  j  ■  L  (>  J-l),  the  neglected 
remainder  product  in  the  denominator  is 


Rem 


Lj»L+l 


i  £  r 


L+l 


( E-2 ) 


This  is  substantially  the  same  as  (E-l),  where  terms  of  the  order 
2L 

of  r  have  been  neglected. 
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